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We derive a first-principles method of determining the van der Waals or Casimir pressure in a 
dissipative and dispersive planar multilayered system by calculating the Maxwell stress tensor in 
a fictitious layer of vacuum, that is eventually made to vanish, introduced in the structure. This 
is illustrated by calculating the van der Waals pressure in a thin film with dissipative properties 
embedded between two semi-infinite media. 



van der Waals (vdW) forces, resulting from the al- 
teration of the quantum and thermal fluctuations of 
the electrodynamic field due to the presence of inter- 
faces, play a significant role in the interactions between 
macroscopic objects at micrometer and nanometer length 
scales. Hamaker was the first to extend the concept of 
London-vdW forces between two atoms to forces between 
macroscopic spheres by pairwise summation of the inter- 
action energy between atoms that constitute the spheres 
PP. Lifshitz, in his seminal work [2], outlined a method 
based on Rytov's theory of fluctuational electrodynam- 
ics [3] for computing the vdW forces between two semi- 
infinite regions separated by a vacuum gap. It required 
the calculation of the average value of the Maxwell stress 
tensor in the vacuum gap. The generalization of Lifshitz's 
method to calculating vdW forces between semi-infinite 
regions separated by dissipative media is not straightfor- 
ward because of the difficulty in defining an electromag- 
netic stress tensor in dissipative media [4 . The goal of 
this paper is to obtain the general theory of vdW pressure 
in arbitrary planar media with dissipative and dispersive 
electromagnetic properties without resorting to defining 
the electromagnetic stress tensor or free energy in any 
material but vacuum. 

An approach proposed by Dzyaloshinskii, Lifshitz, and 
Pitaevskii (DLP from now on) [5], shrouded in the com- 
plicated language of quantum field theory, is the most 
frequently used generalization of Lifshitz's method to cal- 
culate forces between objects separated by absorbing me- 
dia. Even though it has been noted that an expression for 
Maxwell stress tensor for time-varying fields in absorb- 
ing media cannot be expressed in terms of the frequency 
dependent permittivity and permeability alone [4 , DLP 
method effectively reduces to using a "Minkowski-like" 
[6] definition of electromagnetic stress tensor in dissipa- 
tive media. Ninham et al. [7 circumvented the compli- 
cations of the DLP method but, in doing so, had to pos- 
tulate that the free energy of an electromagnetic mode 
at frequency ujj is given by fc#T log [sinh (fvujj /2kBT)\, 
where ks is Boltzmann's constant, 2ith is Planck's con- 
stant, and T is the absolute temperature, even though 
the mode frequencies in dissipative media are, in general, 
complex. It has been argued by Barash and Ginzburg 
[8j [9] that ascribing to each mode a free energy of the 



above-mentioned form is indeed correct. The methods 
of DLP and Barash and Ginzburg are justified on the 
grounds that it is possible to ascribe thermodynamic func- 
tions to electromagnetic fields in equilibrium with matter 
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The relative transparency of the Lifshitz method is ob- 
scured by the complexity of Dzyaloshinskii's formalism or 
by having to define the free energy of each mode, even 
though the final result is a simple generalization of the 
Lifshitz formula. It has been generally regarded that Lif- 
shitz's method, in which the definition of stress tensor 
is above reproach, is incapable of handling dissipative 
media without relying on either of the two generaliza- 
tions [6 . Using the fluctuation-dissipation theorem and 
properties of the dyadic Green's function, we express the 
components of the Maxwell stress tensor in vacuum in 
terms of components of the dyadic Green's function [TO] . 
After a description of a general method to deal with mul- 
tilayered media, we show, using examples of (1) a thin 
film bound by vacuum on both sides, (2) a thin film with 
vacuum on one side and a semi-infinite medium with ar- 
bitrary permittivity and permeability on the other, and 
(3) a thin film bound by semi-infinite media with arbi- 
trary permittivity and permeability, that the expression 
for vdW pressure coincides with that of DLP. 

Let us analyze a general multilayer system, as shown in 
Fig. [T^i, and express the vdW free energy of the system 
in terms of combinations of vdW free energy of smaller 
units. The vdW free energy per unit area of a planar con- 
figuration of N layers (Fig. [IJi) sandwiched between two 
semi-infinite media, medium L to the left and medium 
R to the right, is represented by Ulr (21, ■ • • , zn). Each 
layer is characterized by not only the thickness z^ but also 
the permittivity, and permeability, /!& (both relative 
to that of vacuum). We use the aforementioned notation 
for free energy for its efficiency. If one of the semi-infinite 
media is vacuum, the subscript V is used instead of L or 
R- Ulr (zi, • • • , zn) can be written as a combination of 
three terms: (1) the free energy of the first k layers sand- 
wiched by semi-infinite medium L to the left and vacuum 
to the right of the k th layer, Ulv (zir " ^ z k\ (2) the free 
energy of the remaining N — k layers sandwiched by semi- 
infinite medium R to the right and vacuum to the left of 
the (k + l) th layer, Uvr (^fc+i, * * * , 2/v)> and (3) the work 
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done in bringing the two systems from infinite separation 
to a separation 5 — >• 0. This statement can be written as: 



Ulr (zi, • , z N ) = U LV (zi, • • • ,Zk) + 

5 

U VR (z k+ u • • • ,z N ) + lim / T™ 9 (z v ) dz v , 



(1) 



where T™ 9 (z v ) = T™ 9 (z lr • • • ■ • ,*i>r) is 

the vdW pressure in the vacuum region in Fig. [TJd against 
which work needs to be done to create the N layer sys- 
tem from the two sub-systems. The partial derivative 
dU LR (zi,-'-, z N ) /dz r =p^ (zi,--- , ^at) gives the vdW 



pressure in the r th layer of the N layer system bounded 
by L and i?. For a thin film bounded by two semi-infinite 
regions, we drop the superscript (r) and denote the pres- 
sure simply as Plr- By differentiating Eq. [TJwith respect 
to z r , we obtain the following equation for Pi) R '. 
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One of the first two terms on the rhs of Eq. [2] is zero 
depending on whether l<r<kork + l<r<N. 
Though T™ 9 (z v ) diverges as z~ 3 for z v — >• 0, the quan- 
tity dT^ 9 /dz r is finite as z v V 1 < r < TV, allow- 
ing us to define the partial derivative of the last term in 



Eq. 



as the integral j dT^ 9 /dz r dz v (see supplemen- 



tal information for a justification). T™ 9 (z v ) is obtained 
simply by determining the zz component of the Maxwell 
stress tensor in the vacuum region. Using the proce- 
dure described above, we can write the vdW free energy 
of any N layer medium in terms of Uvv(zi), ^1/(^2), 
• • • , and Uvv(zn), and contributions from terms of the 
5 

form j T™ 9 (z v )dz v , all of which involve calculation of 
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Maxwell stress tensor in vacuum alone. Uyy{z) is noth- 
ing but the vdW free energy to create a thin film of thick- 
ness z in free space. 

We rely on Rytov's theory of fluctua- 
tional electrodynamics to determine the value 
of T™ 9 . The cross-spectral correlations of 

the electric field components can be written 
as (E p (r, uj)E* q (r, 0j))=(2ujfi o G/7r)ImG e pq (r, r) 

(p,q=x,y,z), where 0=(/i<^/2)coth(/iu;/2/c£T), and 
Gp q is the pq component of the electric dyadic 

Green's function [lOl [11]. The spectral correla- 
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tion is defined such that (E p (r,t)E q (r,t)) = fduj 



(E p (r,uj)E*(r,uj)}. Similarly, the cross-spectral cor- 
relations of the magnetic field components are given 
by (ff p (r,o;)fr*(r,a;))=(2a;e o e/7r)ImG^(r,r). The 
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FIG. 1: (a) A multilayer system with N layers between two 
semi-infinite regions L and R. (b) Method of splitting N 
layer mutilayer system into components. The z-axis is per- 
pendicular to the interfaces. 



dyadic Green's functions G e and G m are electro- 
magnetic duals of each other and are solutions of 
\7x\7xG(r,r')-k 2 G(r,r')=I5(r - r'), where J is 
the identity dyad. G e and G m are obtained by en- 
forcing the continuity of: (1) p(r)(nxG e (r,r')), (2) 
nxVxG e (r,r'), (3) s(r)(nxG m (r,r')), and (4) 
nxV 'xG m (r,r') on either side of an interface defined 
by the unit normal vector h at the point r. 

The zz component of the Maxwell stress tensor 
in vacuum can be expressed in terms of G e and 
G m as T zz (r,uj)=(2ujG/7rc 2 )ImG (r,u) QE1 ID], where 
G (r, u>)=Gl z (r, r)-±TrG e (r, r)+G^(r, r)-\TvG™{r, r). 
The average value of the zz component of the Maxwell 
stress tensor, T™ 9 , at any instant of time at position r 
in a vacuum layer is given by: 
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G e and G m are analytic in the upper half plane (UHP) 
by virtue of being response functions. Since G{r,uj) is 
a linear combination of different components of G e and 
G m , it is also analytic in the UHP. We can therefore 
use Lifshitz's technique to replace the integral over uj 
along the real positive frequency axis by a summation 
over Matsubara frequencies on the imaginary frequency 
axis in the UHP as: 



rj-iavg 
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J2'£ G ( r ^n) = k B Tj2' K n (4) 

n=0 n=0 

where, C n =2irnk B T/h, K n =-2&G(r,i£ n )/c 2 , K = 
-lim2^ 2 G(r,^)/c 2 , and ra=0, 1, 2,- • • . The prime (') 

next to ^2 indicates that the n = term is given weight 
0.5. G(r,z£ n ) can be written in terms of the reflection 
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FIG. 2: A four layer system. 



coefficients of plane waves that comprise G e and G m 
[T2] . We now apply this method to calculating the vdW 
pressure in a thin film (indicated by m) bounded by two 
semi-infinite objects, L and R. To do so, we introduce 
a vacuum layer, shown in Fig. |2j in which the Maxwell 
stress tensor will be determined. 

We start with the assertion that the vdW pressure in 
any infinite or semi-infinite planar medium is zero. We 
will show using the following three examples that the 
proposed method is in agreement with the predictions of 
DLP theory for the case of a thin film between two semi- 
infinite objects. K ni from which T™ 9 can be calculated 



using Eq. [4j for the configuration shown in Fig. [2] is 
given by: 
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P=e,i 



where p = e,h refer to the transverse electric and trans- 
verse magnetic polarizations respectively, and 
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k Z v = yJkl + Evil/c 2 , k zm = yjk 2 +£ m £ 2 /c 2 , (6c) 



(6a) 
(6b) 



and similarly for reflection coefficients at other inter- 
faces and wavevectors in other layers. All permit- 
tivities and permeabilities are evaluated at z£ n , n = 
0, 1, 2,- • • . For reflection coefficients, the subscript v will 
be used to denote an interface with vacuum. Since it is 
o^— J oo T^ 9 (z m ^ z v )dz v that will eventually be used in 
calculating vdW pressure, we give below the expression 
fo r f° ^dz 
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Example 1: Vacuum-Thin Film-Vacuum - To find the 
vdW pressure in a thin film of material of thickness £ m , 
we consider a four layer configuration, as shown in Fig. |2j 
with L being replaced with material m, and R being vac- 
uum. If the vdW energy for creating a film of thickness 
z m is Uyy{z m ), the following equation can be written 
for conservation of energy for moving the thin film from 
z v = oo to z v = 5 — » 0: 



Uvv(zi 



o 

+ lim / T^dz v = U , 
s^o J 



(8) 



where, U Q is an arbitrary constant that is the energy 
per unit area of a semi-infinite medium M adjacent to 
a semi-infinite region of vacuum. Differentiation Eq. [8] 
with respect to z m gives the following equation for vdW 
pressure: 
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(9) 



Using Eq. |4j Eq. [7| and Eq. [9j the vdW pressure in a 
thin film of medium m bounded by vacuum is given by 
(see supplemental information for further details): 
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where each integral is evaluated at the Matsubara fre- 
quency £ n = 27rnkBT/h. 

Example 2: Material-Thin Film-Vacuum - To find 
PLv(zm), we consider a four layer configuration, as shown 
in Fig. [2j with R being vacuum. Equation [I] can be 
modified for the four layer system to give the following 
equation for ULv(zm): 

5 

U L v(z m ) = U vv (z m ) + lim / T™ 9 dz v , (11) 



Differentiating Eq. 11 with respect to z m , we obtain 
the following equation for pLv(zm) in terms of pvv(zm). 
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which has been calculated earlier, and T^ 9 : 
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Using the expressions for pvv (Eq. [l0| ) and Eq. [7| we 
obtain the following equation for pLv(zm) (see supple- 
mental information for further details): 
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We can obtain the vdW pressure p^(2; m ) by replacing 



L with i? in Eq. 13 



Example 3: Material-Thin Film-Material - The vdW 
free energy of the system L — m — R is obtained by adding 
to the free energy UyR(z rn ) the work done in moving this 
system from infinite separation to the surface of a semi- 
infinite region of material L. Written as an equation, we 
get: 



U L R(Zm) = U V R(Zr, 



lim / T2Z°dz v , (14) 



along the real frequency axis. The extension to a mul- 
tilayered medium is simply an exercise in determining 
the appropriate reflection and transmission coefficients 

nana. 

We have provided here a transparent formalism for cal- 
culating vdW or Casimir pressure in dissipative and dis- 
persive media that are constituents of planar multilayer 
structures without having to define or calculate the stress 
tensor in such layers. We provide evidence backing the 
generalization of Lifshitz theory of vdW forces without 
relying on quantum field theoretic techniques employed 
by Dzyaloshinskii, Lifshitz, and Piatevskii. These results 
offer further proof of the validity of the Minkowski-like 
stress tensor for calculating vdW forces, at least in planar 
multilayered media. This formalism can be generalized 
to obtain the vdW free energy and pressure of systems 
involving finite sized objects. 

The authors would like to acknowledge correspondence 
and discussions with Prof. Adrian Parsegian, Prof. Rudi 
Podgornik, and Prof. L. Pitaevskii. It was brought to 
our attention that Pitaevskii and Lifshitz had started to 
solve the problem of vdW pressure in dissipative media 
by the method outlined in this paper but put an end to 
it on the suggestion of Landau [T5]! 



and PLR(zm) is given by (see supplemental information 
for further details): 



PLR(Zm) = PVR(Zm) + 
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It can be seen that Eq. 15 for pLR(z m ) is a generalization 
of Eq. [lO] and Eq. [13} Further simplification of Eq. [l5j 
as shown in the supplemental information, results in the 
following expression for Plr' 
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where q = k zm /(€ n y/e m (i€n)/c). Equation [l6| agrees 
with the expression for vdW pressure in a thin film ac- 
cording to DLP [5j [13] . We stress that the method out- 
lined here for calculating vdW pressure is valid irrespec- 
tive of computation of the electromagnetic stress tensor 
by a summation along the imaginary frequency axis or 
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